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Abstract. - Adiabatic quantum charge pumping in graphene nanoribbon double barrier struc- 
tures with armchair and zigzag edges in the resonant transmission regime is analyzed. Using recur- 
sive Green's function method we numerically calculate the pumped charge for pumping contours 
encircling a resonance. We find that for armchair ribbons the whole resonance line contributes 
to the pumping of a single electron (ignoring double spin degeneracy) per cycle through the de- 
vice. The case of zigzag ribbons is more interesting due to zero-conductance resonances. These 
resonances separate the whole resonance line into several parts, each of which corresponds to the 
pumping of a single electron through the device. Moreover, in contrast to armchair ribbons, one 
electron can be pumped from the left lead to the right one or backwards. The current direction 
depends on the particular part of the resonance line encircled by the pumping contour. 
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Introduction. — In recent years graphene has been 
the subject of intense theoretical and experimental re- 
search mainly due to its very peculiar electronic struc- 
ture. Charge carriers in graphene being effectively mass- 
less are well described by the Dirac-like equation in con- 
trast to normal semiconductors with quadratic dispersion 
law. This results in the Klein tunnelling "paradox" pQ, 
half-integer quantum Hall effect [2] observable even at 
room temperature, and other effects [3HS]. Many authors 
consider graphene as a good candidate for spintronics and 
, for future replacement of silicon-based electronics. How- 
ever, gapless nature of an infinite graphene sheet is the ori- 
gin of low on-off current ratios of graphene-based FETs. A 
graphene sheet can be cut to form graphene nanoribbons 
(GNR) with different orientations of edges relative to the 
graphene crystal structure. This leads to the opening of 
finite band gap due to additional transverse confinement 
of the carriers. The band gap is nonzero for both armchair 
(AGNR) and zigzag (ZGNR) ribbons, its value being de- 
pendent on the ribbon's edge type and its width. The 
finite band gap substantially increases on-off current ra- 
tio of GNR-based FETs [SHE] . Experimental prototypes of 
graphene-based FETs were demonstrated to operate under 
tens of GHz frequencies [3[9] and are expected to outper- 
form their silicon-based counterparts. 



A response of a mesoscopic system to a time-dependent 
perturbation has attracted a lot of interest. If two (or 
more) independent parameters (e.g., gate voltages) of a 
mesoscopic system arc adiabatically modulated in time, fi- 
nite dc current through the device can be generated. This 
phenomenon is known as adiabatic quantum pump effect. 
A quantum pump can be used as a quantum standard 
for current if the charge pumped through the system per 
cycle of modulation is quantized |10j . Such quantization 
is naturally achieved in devices operating in the Coulomb 
blockade (CB) regime [11]. But as it turns out, CB is not a 
necessary condition for the charge quantization. Y. Levin- 
son et al. considered |12j a quantum dot separated from 
the leads by two potential barriers, whose heights serve 
as pumping parameters. They argued (in neglecting CB 
effects) that the charge pumped through the dot per cycle 
of modulation is close to a single electron charge when the 
pumping contour encircles the peak of resonance transmis- 
sion [P2lll3| . Clear physical picture of charge loading and 
unloading, explaining this quantization, was later elabo- 
rated in ref. [14]. 

Quantum pump effect in graphene was studied by 
E. Prada et al. using Dirac approximation [T5]. They ar- 
gued that the Klein tunnelling effect has a great impact on 
the properties of graphene-based pumping devices due to 
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Fig. 1: Schematic of the quantum pump ZGNR-based device. 
On-site potentials U\ and XJ% have width of N\ unit cells. The 
right barrier is shifted by half a unit cell so that the whole 
structure is left-right symmetric. The central region comprises 
No unit cells. AGNR case is similar and hence not shown. 
Pumping is achieved by periodical variation of potentials Ui 
and U2. Typical pumping contour is shown on the right. 

the unusual (in comparison with normal devices) contri- 
bution of evanescent modes. Pumping with two potential 
barriers, separated by finite unbiased central region, was 
considered in ref. [16 . It was demonstrated that due to 
the high anisotropy of transmission through a potential 
barrier in graphene both directions of pumping can be re- 
alized for a fixed pumping contour in contrast to normal 
devices. Pumping with a series of barriers was considered 
in a recent paper by Z. Wu et al. [T7] . 

In this paper we study adiabatic quantum pumping in 
graphene nanoribbons with zigzag and armchair edges. 
The main attention is paid to resonant tunnelling regimes 
and the quantization of the pumped charge. Although 
pumping in AGNR device resembles that in conventional 
quantum dots, the ZGNR case is very much different due 
to peculiar tunnelling properties inherent to ZGNR. 

Model. — The geometry of the setup is depicted in 
fig. [TJ The Hamiltonian of the device can be expressed 
using orthogonal nearest-neighbour tight-binding approx- 
imation with one 7r-electron per atom in the form 
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Here (ij) denotes the summation over nearest neighbours 
with hopping parameter t = 2.7 eV. Pumping is achieved 
by periodic variation of two external gate voltages which 
are modelled by adding on-site energies XJ\ and t/jj to the 
diagonal terms of the Hamiltonian. Electron-electron in- 
teractions are neglected and the spin degeneracy factor 
of 2 is omitted for clarity throughout this paper. 

In the adiabatic approximation the charge entering the 
lead a during one pumping cycle with period T is given 
by Brouwer's formula [TH] 
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where e is the electron charge and the emissivity dN a /dUi 



is defined by the expression 
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where Sik is the scattering matrix of the device, and the 
summation over i is restricted to the open channels in the 
lead a. The emissivity dN a /dU is the local partial density 
of states (integrated over the region where potential U is 
applied) associated with the carriers entering the lead a 
regardless of the lead from which they were injected [191 

H]]. 

Expression ([2]) can also be rewritten as the surface in- 
tegral over the pumping contour area A 

Q a = - [ dU 1 dU 2 U a (U u U 2 ), (4) 
n J a 

k iEa 

In the rest of the paper we consider only Q = Qr, and 
Ql = —Qr due to the charge conservation. 

Emissivity can be calculated directly by numerical dif- 
ferentiation of the scattering matrix in eq. ([3]) or using 
the following expression which can be derived within the 
Green's function formalism 1211 



dN a 
~dU~. 
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Here G 7 '^ is the retardered (advanced) Green's function 
of the device and r Q is the line-width function of the 
lead a. The potential profile is described by the diagonal 
matrix Aj with elements (Aj)jj = 1 if the site j belongs to 
the region where Ui is applied, and (A-i)jj — otherwise. 
The Green's function of the device is determined by 

G r = [E - H - Si - X R }-\ (7) 

where Hq is the Hamiltonian of the device region, and 
self-energies and account for semi-infinite left and 
right leads (modelled by graphene ribbons), respectively. 
The Green's function can be calculated either using direct 
matrix inversion (for small devices) or using recursive al- 
gorithms. Self-energies Y<l,r of the leads can be obtained 
using iterative scheme or, for instance, using eigendecom- 
position method [52]. Then the line- width functions are 
given by T a = i(E a — EJ,). Once the Green's function of 
the device is known, the scattering matrix can be obtained 
using the Fisher-Lee relation [23] . 

In this paper we focus our attention on resonant tun- 
nelling regimes with only one open channel in the leads. 
Then expression ([5]) for U(Ui,U 2 ) = Hr(Ui,U 2 ) reduces 
to 

dt* dt dr* dr 



H(Ui,U 2 ) =Im 



(8) 



dUidU 2 dU x dU 2 

where t (r) is the transmission (reflection) amplitude for 
right-moving carriers. We also used the fact that the 5- 
matrix is symmetric in the absence of a magnetic field 

EH. 
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Fig. 2: (Color online) Contour plot of IT(t/i, U 2 ) for 10-AGNR. 
Two pumping contours are shown: AB (red) and CD (purple). 
At the points A and C (B and D) the quasi-bound level in- 
side the device moves below (above) the Fermi level Ep and 
an electron tunnels from the left lead into the central region 
(from the central region into the right lead). Inset: contour 
plot (logarithm scale) of the conductance G{U\,U2) /Go. The 
parameters are: E F = 0.637 eV, iV = 30, N\ = 10. 

From the scattering matrix the conductance G(Ui,U2) 
of the device can be calculated via the Landauer-Biittiker 
formalism: 

G(U 1 ,U2) =G \t\ 2 , G = ^. (9) 

h 

We now proceed with the numerical results. 

Armchair ribbons. — In fig.[5]we plot H(Ui, [7 2 ) and 
G(Ui, U2) for 10-AGNR device (Fermi energy in the leads 
is fixed). The resonance line of the conductance G(U\, U2) 
(see the inset) on the (Ui, U2) -plane corresponds to the 
alignment of the Fermi energy Ep in the leads and the 
energy of a quasi-bound state in the central region, which 
depends on U\ and 1/2- The resonance peak of conduc- 
tance is accompanied by the resonance peak of IL(Ui, J/ 2 )- 
If the pumping contour encloses this peak, the pumped 
charge is quantized [T2HT4] . We stress that this quan- 
tization has no relation to the CB as far as we neglect 
electron-electron interactions. The quantization is only 
approximate because n({7i, U2) remains finite outside any 
finite pumping contour. 

Charge quantization can be intuitively explained as the 
loading and unloading of one electron into/out of the de- 
vice (see |14j for details). Consider the pumping contour 
AB shown in fig. [2] It has two resonance points A and 
B, where G(U\ 1 U2) attains its maximum value (but still 
G(Ui,U~2) <C Go at these points). When the point A is 
crossed, the quasi-bound level inside the device moves be- 




Fig. 3: (Color online) The pumped charge per cycle Q for 10- 
AGNR as a function of the Fermi energy Ep in the leads for 
fixed pumping contours: AB (red solid line) and CD (purple 
dashed line). The transmission maxima (with fixed Ui — U2 — 
0.36 eV) are shown with the black arrows. The parameters of 
the device are the same as in fig. [2] 

low the Fermi level in the leads and an electron tunnels 
from the left lead into the central region. The tunnelling 
probability from the right lead is much smaller because 
the right barrier is much higher. At the point B the sit- 
uation is reversed: the quasi-bound level moves up and 
an electron tunnels into the right lead. Hence, a single 
electron is transported from the left to the right per one 
pumping cycle, and therefore Q = +e. 

Let us now fix two pumping contours (AB and CD) 
and consider the dependence of Q on the Fermi energy 
Ep. If Ep is varied, a resonance line on the (Ui, C/2)-plane 
appears whenever Ep becomes equal to the energy of some 
quasi-bound state. When the Fermi level is far from any 
such state, no resonance is observed and G(Ui, U2) <C Go 
everywhere inside the contour AB. One can see from fig. [3] 
that for some values of Ep the pumped charge for both 
contours tends to the value +e. This happens each time 
when the pumping contour encloses the resonance peak of 
1/2)- One can expect less accuracy of quantization 
for the contour CD than for the AB one because the former 
is much smaller and it encloses lesser part of the resonance 
peak. Due to the same reason the peaks corresponding 
to the contour CD are narrower. Each peak of Q(Ep) 
is accompanied by a peak of the transmission (for fixed 
potential barriers) shown in the same figure. 

It is known that ./V-AGNR becomes metallic whenever 
N = 3n — 1 for some integer n. The lowest mode has 
zero transverse momentum and experiences perfect Klein 
tunnelling through potential barriers and, hence, cannot 
be pumped |15) . However, the finite gap opens upon in- 
troduction of nonzero next-nearest-neighbour hopping pa- 
rameter into the simple nearest-neighbour tight-binding 
model (fTJ) [25,26 . The results are then qualitatively sim- 
ilar to that presented above. 

Zigzag ribbons. — The situation is qualitatively dif- 
ferent for ZGNR. We plot H(U U U 2 ) and G{U U U 2 ) in 
fig. |U One can observe that the resonance line is sepa- 
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Fig. 4: (Color online) Contour plot of U(Ui, U 2 ) for 10-ZGNR. 
Three pumping contours are shown: BC (red), AB (blue) and 
EF (purple). BC (AB) encloses a positive (negative) peak and 
pumps Q = +e ( — e). Inset: contour plot (logarithmic scale) 
of the conductance G(TJ\, U2)/Gq. At the points A and C (B 
and D) the left (right) barrier is opaque. The parameters are: 
E F = 0.36 eV, Ao = 30, Ni = 8. 



rated into several parts in which H(Ui,U2) has either a 
positive or a negative peak. If a pumping contour encir- 
cles a positive (negative) peak, the charge Q = +e (— e) 
is pumped per cycle. When we follow along the resonance 
line, the peak sign of H(Ui, U2) alternates. This peculiar 
behaviour can be understood by considering the inset in 
fig. E] where the contour plot of conductance G(U\, U2) is 
shown. As can be seen from the figure, the conductance 
vanishes on certain lines on the (Ui , U2) plane. These zero- 
conductance resonances (dips) [27 [29] are associated with 
the formation of discrete quantum levels in the barrier re- 
gion. Their existence in the case of ZGNR with single 
applied potential barrier was previously demonstrated nu- 
merically in ref. |29j . The zero-conductance lines can be 
approximated by straight lines (U1.2 = const). Vertical 
(horizontal) lines correspond to the dips associated with 
left (right) barrier, and the barriers act effectively inde- 
pendently. In close vicinity of the points of intersection 
of vertical and horizontal lines the approximation of in- 
dependent barriers breaks down. Each zero-conductance 
line is accompanied by a 7r-phase jump of the transmission 
amplitude [25] . 

Now let us follow the contour BC. The point B (C) be- 
longs to the horizontal (vertical) zero-conductance line. 
When we cross the point B, the height of the left barrier 
U\ decreases and so the quasi-bound level inside the de- 
vice moves below the Fermi level in the leads. The right 
barrier is opaque and an electron tunnels from the left 
lead into the central region. At the point C the height 
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Fig. 5: (Color online) The pumped charge per cycle Q for 
10-ZGNR as a function of the Fermi energy Ef in the leads for 
fixed pumping contours: BC (red solid line), EF (purple dashed 
line) and AB (blue dotted line). The pumped charge for 9- 
ZGNR for the contour BC is also shown with the black dashed- 
dotted line. The transmission maxima (with fixed Ui = U2 = 
0.45 eV) are shown with the black arrows. The parameters of 
the device are the same as in fig. [4] 

of the right barrier U2 increases and the level moves up 
forcing an electron to tunnel into the right lead because 
the left barrier is opaque now. Hence, this pumping cycle 
transports a single electron from the left to the right, and 
Q = +e. The situation is reminiscent of that for AGNR. 

Finially we analyze the contour AB. At the point A (this 
point lies on the vertical line) the level moves down. But 
now the left barrier is opaque and an electron tunnels from 
the right lead. At the point B an electron tunnels into the 
left lead. Hence, the situation is reversed with respect to 
the contour AB: an electron is transported from the right 
to the left and thus Q = — e. The same argument holds 
for the contour CD. 

In fig. [5] we plot the dependence of the pumped charge Q 
on the Fermi energy Ep for the pumping contours AB, BC 
and EF. For Ef = 0.36 eV Qab approaches the value — e 
and Qbc approaches the value +e in agreement with the 
above considerations. It is interesting to note that the 
contour BC pumps the charge — e for Ep — 0.28 eV. This 
can be understood by a similar argument. The relative 
positions of the pumping contour, the resonance line and 
zero-conductance lines is such that the point B now lies 
on the vertical zero-conductance line and the point C lies 
on the horizontal one. Thus, the picture is qualitatively 
similar to that for the contour AB in fig.|U The left barrier 
is opaque at the point B and the right one is opaque at the 
point C. The reasoning proceeds as for the contour AB, 
resulting in Q — —e. Just like in the AGNR case, peaks 
of Q{Ep) and peaks of the transmission come in pairs. 

For anti-zigzag ribbons (7V-ZGNR with odd N) with the 
selected form of the pumping potentials the transmission 
through the device does not exhibit resonant tunnelling 
behaviour. Hence, we do not expect the charge quantiza- 
tion in this case. This is confirmed in fig. [SJ where the 
pumped charge for 9-ZGNR is plotted. 
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Fig. 6: (Color online) The pumped charge per cycle Q for 
12-ZGNR with smoothed pumping potentials as a function of 
the Fermi energy Ef in the leads for fixed pumping contours: 
BC (red solid line) and AB (blue dotted line). Inset: schematic 
of pumping potentials profile. The parameters are: iVo = 38, 
Ni = 7,a = 0.2, AB: Ui = 0.4 H- 0.43 eV, U 2 = 0.52-^0.64 eV, 
BC: Ui = 0.43 "J - 0.53 eV, Ui = 0.43 ~ 0.53 eV. 

Finally, we consider the approximation of rectangular 
potential barriers. The profile of pumping potentials gen- 
erated by metallic gates should be determined via self- 
consistent solution of Schrodinger and Poisson equations. 
However, it was shown numerically in ref. [29 that the 
transmission through a smooth potential barrier in ZGNR 
also demonstrates zero-conductance resonances. Hence, it 
is reasonable to expect the same qualitative pumping be- 
haviour for more realistic smoothed potentials. To demon- 
strate this, we employ smoothed potential barriers 

( U a e~ a i\ i<0 
U(j) = \ U , < j < JVi (10) 

I u e- a ^- N ^ 2 , j>m, 

where j is the coordinate (in the unit cells) along the rib- 
bon, Uo is the amplitude of a pumping potential and a 
is the smoothing parameter. The numerical results pre- 
sented in fig. [HI are in agreement with the expectations. 

Conclusion. — We have analyzed adiabatic quantum 
charge pumping of nonintcracting electrons in AGNR and 
ZGNR double barrier structures in the resonant transmis- 
sion regime. 

We consider AGNR ribbons first. In this case the whole 
resonance line contributes to the pumping of a single elec- 
tron (ignoring double spin degeneracy) per cycle through 
the device. This picture is reminiscent of that of a sim- 
ple ID double barrier structure (a quantum dot separated 
from the leads by two point contacts with variable conduc- 
tances). When the pumping contour encircles a large part 
of the resonance line, the current is quantized. The direc- 
tion of current is defined by the direction of the pumping 
cycle. 

The existence of zero-conductance resonances makes the 
ZGNR case more complicated and qualitatively different. 
These resonances separate the whole resonance line into 
several parts, each of which corresponds to the pumping 



of a single electron through the device. An electron can 
be pumped from the left lead to the right one or back- 
wards. The current direction depends on the particular 
part of the resonance line encircled by the pumping con- 
tour. This behaviour stems from the zero-conductance 
resonances inherent to locally gated ZGNRs. 

Two points toward the experimental realization of 
pumping devices described above have to be mentioned. 
First, the size of the devices studied in this paper was pri- 
marily limited by the computational cost, but the same 
picture should be observable for longer and wider devices 
provided that the quasi-bound states in the central region 
are well resolved. Second, the edge shape is expected to be 
crucial for the observation of peculiar pumping behaviour 
of ZGNR-based devices. Although the precise control over 
the edge type still remains a challenging task, the technol- 
ogy advances very fast [30H32] and the controlled pattern- 
ing techniques of GNRs with various sizes and edge types 
might become available in the near future. 
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